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Abstract 

In this paper, to obtain an analogy between the curved spaces and the hnear optics, 
we expand the idea of Ref.[Tl [5] to the multilayer films. We investigate effects of 
thickness and index of refraction of the films on the Lorentzian transformations. In 
addition, by using the multilayer films, we suggest very simple experimental set-up 
which can serve as an analogue computer for testing special relativity. Finally, we 
draw an analogy between the Rindler space, as an example of the curved spaces, and 
a suitable multilayer film. 
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1 Introduction 

In recent years, analogies have taken on a major role in physics and mathematics [3ll31[5]. 
They provide the new ways of looking at problems that permit cross-fertilization of ideas 
among different branches of science. Analogy methods have practical worth from two point 
of views: First, they provide new light on perplexing theoretical questions OH]. Second, 
some of them are interesting for experimental set-up [3]. 

In this contribution we use the method developed in [B [2]) in a multilayer stack to 
show how we can set up an experimental model for special relativity and Rindler space as 
a special case of general relativity. 

In the present section we review some of the relevant results achieved in p4 i2j . As it 
is well known, SL{2, C) group is locally isomorphic (or a 2 to 1 homomorphic) to Lorentz 
group, SO{l,3). In fact, the transfer matrix, used in the classical optics for description of 
the interaction of light and matter, can be considered as an element of SL{2, C) (see [6] 
and references therein). 
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In ref. [2], authors gave a geometrical interpretation for the transfer matrix of an 
absorbing system. In fact, they considered an isolated absorbing film embedded between 
two identical ambient I and substrate // media. A monochromatic, linear polarized plane 
wave falls from the ambient, as well as a plane wave of the same frequency and polarization 
falls from the substrate. If the amplitudes of the field are characterized by a column vector, 

HI) 

where E_ is denoted backward (forward)-travelling plane wave, and the amplitudes 

of the ambient and substrate are related by 

Ej = M Ejj, (2) 

where M is the transfer matrix and is defined as, 

/ T'^-R^ R \ 
M = 3j T , (3) 

\ T T / 

and T and R are the transmission and reflection coefficients, respectively. 
By defining the coherency matrix [7] as 

^=[ eJx \eJ j' 

for both / and //, it is easy to show that, 

8i = M8nM\ (5) 
It is useful to identify the variables x^, 

= ^rr(£:a'^), (6) 

as the coordinates of a (1 + 3)-dimensional Minkowski space, in which = 1 (the identity) 
and (T*'s are the standard Pauli matrices. The conjugation ([5]) induces a transformation on 
the variables of the form 

= Tji{M) x'ij, (7) 

where, 

r^(M) = ^Tria'^Ma^M^). (8) 

It is easy to see that the transformation T(M) is a Lorentz transformation [8j which is 
generated by M and —M. Thus it is a 2 — )• 1 homomorphism map between SL(2, C) and 
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50(1,3). 



We use the model which is developed in [U [2] , and consider a monochromatic, linearly 
polarized plane electromagnetic wave incident on a non-absorption multilayer thin film. In 
this model basis coordinates are defined by the components of the electromagnetic wave, 
i.e., E and H. 

The organization of this paper is as follows. In section [2] we review the transfer matrix 
method for description of the multilayer films. In section [3] we introduce spinor map 
between 5L(2,C) and S'0(l,3), and we draw analogy between an EM field which passes 
through a multilayer film and transformation of coordinates under the Lorentzian group. 
In addition, two special multilayer, that is, quarter-wave films and a multilayer 

with very low thickness films, are studied in detail. In section HJ by using the multilayer 
films, we suggest an experimental set-up for investigation of two kinematical properties, 
that is, time delay and length contraction as some of predictions of the special relativity. 
In section [5] we introduce an analogy between the Rindler space, as a curved space, and a 
suitable multilayer film. Finally, a summary and conclusion remarks are given in section 

El 



2 Matrix theory for description of the Multilayer films 

We consider an electromagnetic radiation falling onto a structure consisting of a thin 
film with uniform thickness t and index of refraction n. This isolated (homogeneous and 
isotropic) film, is embedded between two media, / and //, with indexes of refraction, n/ 
and njj, respectively. We consider a monochromatic, linearly polarized plane falling wave 
which is injected into the thin film from the above with amplitude Ej and Hj. The field 
after passing the end layer is characterized by Ejj and Hjj, as it is shown in Figure ([T|). We 
use transfer matrix method which is a method used in optics to analyse the propagation 
of electromagnetic through a stratified (layered) medium [71 [9] . 

The transfer-matrix method is based on the fact that, according to Maxwell's equa- 
tions, there are simple continuity conditions for tangential component of electric field and 
normal component of magnetic field across boundaries. 

In addition, this method can be developed to the multilayer films. [9j. For a stack with 
m layers, the transfer-matrix can be defined as a product of matrices given by 

M = MmMm-i ■ ■ ■ Ml, (9) 

where the transfer- matrix of the j'*'* layer, Mj, is given by 

(r . sin Sj \ 

, (10) 
ir]j sin 5j cos 6j I 
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where Sj = ^tjuj cos{9j) and r]j for s-polarization and p-polarization is characterized by 
rjj = njCos{6j) and rjj = nj/cos{6j) respectively [9j. It is obvious that the angle 9j is 
related to the angle of incidence by Snell's law, that is, 



rij sm Oj = no sm Oq. 



(11) 



3 Lorentzian interpretation of the transfer matrix 

In this section, we use the spinor map, from 5L(2,C), onto the Lorentz group, S0{1,3). 
We define a coherence matrix [7] as follows: 



£ 



E*H 

EH* 



(12) 



for both / and //, so that £j and Sji are related acording to the Eq. ([5]). We explicitly 
define the coordinates x'^, by using equation 



X 



■° = i(|S|2 + |F|2), x^ = ReiE*H), 



X 



lm{E*H), 



- \H\'^) 



(13) 



It is worth to note that these coordinates behave as a light-like interval: 



(x0)2 - {x'j' - {x'Y - {x-') 



„1\2 



2\2 



0, 



(14) 



Therefore, the conjugation relation ^ induces a transformation on the variables x^ in 
the form of equation ^ . 

However, we can decompose the transfer matrix of a thin film, according to equations 
Q and (do]) with j = 1, as 



Vv 







(15) 



Thus, under the homomorphism ([8]), B generates a boost of velocity P = v {'m natural 
units, where c = 1), while TZ induces a pure spatial rotation. 

To gain physical insights, we write explicitly the equation ([8]) for B and TZ. Under the 
homomorphism TZ generates the following matrix 



T{TZ) = R{25) 



/ 1 \ 

10 

cos 25 sin 2(5 

\ -sin2<5 cos 2(5 / 



(16) 
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which is a rotation in the plane — of angle twice the phase of the transmission 
coefficient. 

Furthermore, diagonal matrix B, under the homomorphism ([8]), generates a boost in 
the plane — [lOj , 



r(e) = L(o,o,/3) 



/ coshC sinh^ \ 

10 

1 

y sinhC cosh(" ) 



(17) 



where rapidity C is obtained from the relations 

coshC = J(r/ + -), (18) 
sinhC= J(r?-i), (19) 

and velocity (3 is defined by 

iv--) 

/3 = tanhC = — (20) 

+ 

In Figure ([2]) , we have plotted the variation of f3 with respect to the 6 for different values 
of the index of refraction, n, when the electromagnetic field have s-, and p-polarization. 
It is seen that for every values of the index of refraction, /3 has maximum or minimum in 
6 = 0. In addition, for a fixed 9, the /3 is increased by increasing the refraction index, for 
both S-, and p-polarization. However, for the large values of the refraction index, n ^ 1, 
for both polarization, the result is the same, as shown in Figure ([2]). Finally, by using 
equations p^ . (fT6|) and p!7|) we obtain T(M) as 

T{M) = L(0, 0, -P)R{2u;t)L{0, 0, /3) (21) 

where uj = 6/t. Also, with some calculations, we can rewrite equation (f2T]) . in the following 
form [Tn[T2]. 

T{M) = R{2ujt)L{0, P sm{2ut), -/3 cos{2ujt))L{0,0, P). (22) 
In addition, by using Einstein rule for the composition of velocities |13j . 

/3i e /32 = ^ |/3i + 02 + (/3i X (/3i X /^a)) I , (23) 

1 + ^ + 7/3i J 

we can write T{M) as a rotation and a boost, 

r(M) = i?(2a;t)L(0,/3^,/3^), (24) 
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where 



, _ /3sin(2a;t) ^'^ cos{2ujt) 



= —T— ^(cos(2a;t +- 7==]. 25 

'^^ l-/32cos(2a;t) V ^ ' i + ^/T^J ^ ' 

In Figures ([3]) and ([!]), we have plotted the variations of /32, and = \/ P'2 + /^a^ 
with respect to t for n = 2.00049 and 9 = it/9, for s-polarization and p-polarization, 
respectively. As it is seen, /32, (3'^ and depend periodically on t. 

We now turn to the multi-layer stack with ni films. By employing the method which is 
used to obtain equation ([22]) . and by using equation (p^ . we can write T{Ai) as 

T{M) = R{2nt)L{0,B^^''\B^3^^) ■ ■ ■ L{0,B^^\b^^''), (26) 

where = YlT^i "^^^^ 

iS^^) = /3^^)cos(2f]WtW) -/3j)sin(2f](^)t(^)), (27) 

i3f = /3j)sin(20«t«) +/35^^cos(20«t»), (28) 

with O(^) = and 

(i) ^ /j-sin(2a;-t-) (/gQ^ cos(2u;-t') 

^2 1 - (/3^)2 cos{2uH^) ^ 1 + yr^G^^' ^ ' 

- /^^^^ ( ,,,r,J^U^). I (/3(-))^sin2(2^(-)tW) x 

^3 - i_(/3«)2cos(2a;(%«)r°'^^'^ * ^+ i + ^ i_(/3W)2 J' ^^"^ 



Also, by using Einstein addition rule (I23p . we can rewrite equation (126p as 

T(X) = R{2nt)L{0,Bi^\Bi'^'^), (31) 
where ;B2"^^ and B^^^ are given by applying the rule (j23p consecutively. 



3.1 A multilayer with quarter-wave films 

In this section, we consider a multilayer with quarter- wave films [14] . that is tjUj = A/4 
and assume that 9q = 0, so that cos 6j = and sinJj = 1. As an example, for a stack 
consists of two layers, matrix M, is given by, 

/ \ 

M = _^ , (32) 
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which under the homomorphism ([8]), generates a boost as [in 



r(M)^L(0,0,/3fJ 



/ coshcS^ sinhcS^ \ 

10 

1 

y sinhCg^^ coshCqS / 



(2) 

where velocity /3qw is given by 



/3g=tanhC£) 



^ 712 \4 



qw 



l + (^^)4' 



(33) 



(34) 



(2) 

and rapidity is obtained from the relations 

1 //'^l ^ 



coshe' = i((=i)'+(=^)=) 

^ Z\ no J 



n2 ni 

-nhC<S = ^((^)^-(22f 

^ 2\n2 ni 



(35) 
(36) 



In general, for a stack consists of 2m quarter-wave layers, in which the thickness of j 



th 



layer is tj = j = 1, 2, • • • , 2m, after some simple calculation, we can achieve T{M2m) 



as. 



/ coshcj^'") sinhcS"*) \ 



r(M2„) = L(0,0,/3(^'")) 








1 
1 








y sinhCgS"^ coshCg^™^ 



(37) 



where velocity Pqw^^ is calculated by the rule (f23]l . As an example for 2m = 4, we have 



qw 



(38) 



3.2 A multilayer with very low thickness 

In this section we study a multilayer, with m layers, in such a way that t^^^^'s are near to 
the unit, ~ 1. Thus, by using equation (j20p and keeping only first-order term of r/^*), 
we obtain = — 1 + 77^*^ and therefore the Lorentzian transformations are replaced by 
the Galilean transformations. Thus, by these assumption, we can rewrite the equation (I3ip 
as the following relation: 



T{M) = R{2nt)L{0,Bi^\Bi'^''] 



(39) 
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where 

m 

4"^) = ^/3«sin[2(^«+17«)t«], 

i=l 
m 

Bf^^ = J];/3«cos[2(a;(^) (40) 

i=l 

In addition, when we assume t is very thin, from equation (j2ip . by keeping only first- 
order terms in t, the rotation R{2ujt) appears as a pure boost, -L(0, 0,/3), which is followed 
by another pure boost in another direction, L(0, 2a (it, 0), so that the result is like a Thomas 
precession. That is, we can write |101 115]. 

R{2ujdt) = L(0, -2adt, -I3)L{0, 2adt, 0)L(0, 0, /3), (41) 

where u) = ^ncosOu and a = — — is the acceleration along the second axis, x^. 



4 Experimental Set-up of Time Delay and Length Contrac- 
tion by Thin Films 

We can now propose some specific experiments bases on the properties discussed in the 
preceding sections. The basic experimental set-up consists essentially of a monochromatic 
plane wave, linear polarized {s or p) with the specific values of \E\ and \H\. 

We use a stack consists of 2 quarter- wave layers. According to the equation ([33]) . 
the coordinates and under the action of transformation matrix remain invariant. 
Therefore, it is possible to consider T{M) as a boost which belongs to S0{1, 1): 



To show how the experimental set-up works, let us suppose that the electromagnetic 
fields, Ein and Hin, fall onto the multi-layer. From equation (jl3p we have 

Q 1 2 2 3 ^2 2 

and 

= k\Eout\' + \Hout\') = k^fmn\' + {^mn\'}, 

2 2 n2 ni 

4ut = \{\Eout\'-\HU) = \{-)\\Ein?-C^f\mn\']. (44) 

2 2 n2 ni 

These equations clearly show that there are many possible experimental methods to choose 
the input fields in such a way that two different points in a reference frame may be appear 
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in the other frame as either with the same temporal or same spatial coordinate. 

We consider two experimental set-ups with different inputs -E'inQ) and -f^m(j) aiid cor- 
responding to outputs -E'oMt(j) and H^^^f^j-^ ,where j = 1,2. Thus, it is easy to show that 
^L(i) = ^L(2) wlien 

\Ein(2)?' - |-E'm(l)P = ( — )^{|-f^m(l)P " |-f^m(2)P}i (45) 

and = x3^,(2) ^^en 

\Ein(2)\^ - |-E'm(l)l^ = ( — )'^{|-f^m{2)|^ " |-f^m(l)P}- (46) 

rii 

By using equation the first set-up can be used to demonstrate the length contrac- 
tion, 

lout = x1ut{2) ~ ^o«t(l) = ~(^L(2) ~ ^L(l))- (47) 

In order to put forward the phenomenon of the time dilation, let us consider another 
experimental set-up in which Hin(i), -£'m(2) and -ffj„(2) satisfy the equation (f46l) . 

Therefore, measurement of Atout = x^out{2) ~ ^oMt(i) '^^ regarded as a time interval 
measurement. Therefore, we obtain 

^tout = X^out{2) ~ ^OMt(l) = ~(^L(2) ~ (48) 

In summary, we have proposed a series of specific experiments by using a multilayer, 
which contains two films and the thickness of their films are given by tjrij = A/4, that 
clearly demonstrate the most important predictions of the Lorentz kinematics. 



5 A Multi-layer as a Rindler Space 

Maxwell's equations in general coordinates are written as [5] 



i^gg^'E,), = {^g'^B,), = 0, (49) 

£0 



and 



dj^g'^B,) IdU/gg^ 
[ijk] Ekj = — ^^-^^ , [ijk] Hkj = — ^-^^ + no^/gJ , (50) 

where permutation symbol, [ijk], is defined by 

{+1, if ijk is an even permutation of 123, 
— 1, if ijk is an odd permutation of 123, (51) 
0, otherwise. 
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and the Einstein summation convention is used. As it is well known in Optics, the Maxwell's 
equations in empty-curved space can be replaced by Maxwell's equations in the flat space 
which is filled by dielectric media with the following parameters [5], 

e'^ =^i'= = ±^g'K (52) 

Therefore, we can see a correspondence between geometry and matter, as in the general 
relativity. Thus, in this case, the index of refraction is obtained by [5j 



— I £23^31 — £2ie33 £11^33 — £l3e31 £13^21 — £13^23 | (53) 
\ £21^12 — £222 

Now, if we consider the metric, 

ds^ = -dt^ + N^{x){dx^ + dy^ + dz^), (54) 

by using equations (j52]) and ([53]) . the index of refraction is given by n{x) = [N{x)]~^. 

On the other hand, the principle of equivalence make the gravitational field connect 
to the acceleration. Therefore, to study gravitational effects, we can generalize concept 
of Lorentz transformations so as to include observers who are moving with a constant 
acceleration with respect to an inertial frame. This metric is called Rindler metric and 
described as [16], 

ds^ = -{g xfdt^ + dx^ + dy^ + dz^. (55) 




It is easy to see that the metrics (|54p and (|55p are conformal. Accordingly, by using the 
equation (153p , we obtain the refractive index of the analogue media with the Rindler space 
as the following relation: 

n = g X. (56) 

Thus, for analogue Rindler space we assume a multilayer, for which the refractive index of 
the j^^ layer is characterized by, 

nj = g{j - l)t. (57) 

However, up to the first order approximation in /J*^*) and by using equation ()57p . for s- and 
p-polarization, when = 0, we can rewrite equation ([5U|) as, 

T{M) = L{0,M2,B3) (58) 

where 

_ ^ (3(j_-i)tl-± ^ ,2r „,2 
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and 



Also, we can consider a stack consists of the quarter-wavelength layers, i.e., tjuj = A/4. 
Therefore, by using the equation (|56p. the thickness of j'*'' layer is 



*i = + - a;^ = ^ ti. (61) 



After some calculation, we can obtain thickness of the j*'^ layer as. 



i. = y^(x/7-v7^), (62) 
and the index of refraction of this layer is obtained by 

r^j = \l^{y^+VI^)- (63) 

We can now construct an analogy between the Rindler space and a multilayer with 2m 
quarter-wavelength layers, in which the thickness of j*^ layer is given by equation (j62p and 
its refraction index is given by equation (j63p . Accordingly we obtain , 



p _ Vj + 1 + Vj . . 

^Multilyer — ^ _^ ^JJ^ J ^ > 

In Refs. [T71 [H], electrodynamics in the Rindler space had been studied by using a 
phenomenological approach. Accordingly, the component of the electric field, near the 
horizon, x = 0, reads as follow: 

^Rindler(a;) = Eq Ki2n/x{^x), (65) 

where Ki2Tr/x{2TTx/X) is the modified Bessel function of order iln/X. 

Figure ([5]) displays a comparison between the output values of the electric fields are 
given by equation (f64l) and equation (|65]) . i.e.. 



— Abs(|£^|Multilayer " |-E|Rmdler)) 

in terms of x = i / and for different values of A. As it is seen, by increasing x (increasing 



2g-J 

the number of layers), AE approaches to zero. In addition, AE is decreased more rapidly 
by decreasing the A. 
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6 conclusion 



This paper is an example of analogue models of Einstein Relativity. Firstly, we have used 
a geometrical tool to analyze a multilayer film in a concise way that is closely related to 
the special relativity. Then, this method was made an analogy between multilayer stack 
and curved space in the general relativity. In fact, in the first step, we have seen how 
the transfer matrix can be interpreted as multiplication of a boost, a rotation and another 
boost. 

In addition, in this case, we have suggested an experimental set-up which tests the most 
important predictions of the Lorentz kinematics, i.e., time delay and length contraction. 
However, by using this isomorphism, apart from a relativistic presentation of the topic, 
which has interest in its own, we suggested an experimental set-up which is, in theoretical 
view, so simple and understandable, and, in practical level, was so easy to accomplish, in 
contrast to the similar experimental set-up |19j . 

Finally, we have employed analogue method in replacing Maxwell's equations in empty- 
curved space to Maxwell's equations in the flat space where the later is filled by a media. 
Then, we have applied this method to the Rindler frame, which by using principle of 
equivalence, is locally equivalent to gravitational fields. This means that this method 
makes an analogy between accelerated frame in general relativity and a suitable multilayer 
films, with specific values of the refractive indexes and thickness of the films. 
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Figure 1: Scheme of the input and output field in a thin film. 



Figure 2: (Colour online) The variation of /?, for s and p-polarization, with respect to 9 
for different values of the index of refraction, n = 1 are plotted by red and blue line, n = 2 
by red and blue dashing line and n = 10 by red and blue dot line, for s and p-polarization, 
respectively. 




Figure 3: (Colour online) The variation of /Jg by red line, (3'^ by blue dot line, and \/3'\ 
by black dashing line with respect to t for the fixed index of refraction, n = 2.419, for 
s-polarization. 
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Figure 4: (Colour online) The variation of /Jg by red line, /Jg by blue dot line, and 
by black dashing line with respect to t for the fixed index of refraction, n = 2.419, for 
p-polarization 



0.40 - 



\ o.oez \v 



-J with g = 1 and 



Figure 5: (Colour online) The variation of AE with respect to x = 

2m = 60 for A = 0.06199 x 10"^, A = 0.06199 x 10"^ and A = 0.06199 x 10"^ by green ( 
left), red ( middle) and blue (right) points, respectively. 
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